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IMPROVED ESTIMATES OF HF RADIOWAVE FIELD STRENGTH 
NEAR A LAND-SEA INTERFACE 


R. M. Jones 


ABSTRACT 


Results of other researchers are extended to give the detailed 
behavior of mixed-path (land-sea or sea-land) propagation near the 
land-sea interface. It is shown that the mechanism for groundwave- 
mode coupling at the land-sea interface is edge-diffraction, and 
that near the shoreline the direct edge-diffracted wave can make a 
significant contribution to the field. It is further argued that, 
for sea-to-land propagation for an elevated observer, the sea-tvpe 
groundwave modes extend over land, and there is a shadow zone for 
the land-type groundwave modes excited at the shoreline.  Corre- 
sponding results hold for land-to-sea propagation. 


1. INTRODUCTION 


When we measure ocean wave parameters with our Coastal Ocean Dynamics 
Applications Radar (CODAR) (Barrick et al., 1977) the propagation of the radio 
waves from the land-based radar to the ocean target and back depends on how 
radio waves propagate across a land-sea boundary. Of particular interest is 
how the echo strength depends on the location of the transmitting and receiving 


antennas relative to the land-sea boundary. 


Figure 1 shows the geometry of the problem. For large enough distances, 
the propagation is by groundwave propagation rather than line-of-sight rays. 
Because the propagation is reciprocal, I consider here only one-way propagation 
from the scattering patch on the sea surface to the radar on the land. For 
the round-trip path, the amplitude (in decibels) and the phase are simply 
doubled. To increase the generality, however, I consider that both the source 
(over the sea) and the observer (over the land) are elevated. I consider 


special cases later. 


The mixed-path problem in radio wave propagation over the earth has re- 
ceived a great deal of attention. Using the compensation theorem, Wait (1961) 
has calculated coupling coefficients among the groundwave modes at the land- 
sea boundary. Wait and Walters (1963) and Hill and Wait (1981) have calculated 
some numerical examples. Figure 2 shows an example of such calculations for 


conditions relevant to our CODAR measurements. 


Free space 


Source 


Observer 
AE N 
no ÄN 


4 


Sea Land 


£2, U0;0 2 €4; Ho,0 4 


Figure 1. Geometry for propagation over a mixed land-sea path. The source ts 
above the sea (charactertzed by a permittivitu € 9 and conductivity o., or 
a surface impedance z o> and the observer is above land (characterized by a 
permittivity €, and conductivity Gs ora surface tmpedance 2_). The con- 


ductivity of free space is zero, and u, ts the magnetic permeability through 
the medtum. 
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Figure 2. Mixed-path groundwave propagation using the groundwave-mode conver- 
sion formulas of Hill and Wait (1981). The field has been normaltzed to 
the free-space field above a flat, perfect conductor. The source is a 
vertical electric dipole on the surface of the sea 20 km from the shoreline. 
The frequency is 30 MHz. The calculations were made with the program 
MPHGAIN f'rom David Hill (National Bureau of Standards; Boulder, UU Tre 
discontinuity in the field above the shoreline is not considered realtstte. 
The soltd lines are for 200 groundwave modes, the dashed lines for 500 modes. 


One of the most striking features of Fig. 2 is the discontinuitv in the 
field above the land-sea interface. The calculations presented in Fig. 2 use 
200 groundwave modes for both the land and sea. The effect of increasing the 
number of modes to 500 is also shown and is not qualitativelv different. It is 
doubtful that increasing the number of groundwave modes even further would vield 
a continuous solution. The solution shown in Fig. 2 is probably correct for an 
observer far enough from the land-sea interface; it is probably sufficient to 
Smooth the curves by hand to make them continuous (as was done in the figures 
of Hill and Wait, 1981) if one is not especially interested in the exact 
behavior close to the interface. However, for our CODAR application where our 
radar is very close to the shoreline, how one smooths would make a large 


difference on the signal-strength estimate. 


This report more accurately estimates the variation in signal strength 


near the land-sea interface. 


2. SUMMARY OF RESULTS 


i3 Edge-diffraction is the mechanism for coupling of groundwave modes at a 
land-sea boundary. Combining the complex-height, creeping-ray representa- 
tion of groundwave modes with an edge-diffraction coefficient of Kaminetzky 
and Keller (1972) gives the mixed-path groundwave propagation formulas 


of Hill and Wait (1981). 


2) For propagation from sea to land, the sea-type groundwave modes do not Stop 
at the shoreline, but continue over land for some distance for an elevated 
observer. The higher the observer, the farther the modes extend over the 
land. The lower order modes continue farther than higher order modes 
before they are cut off. For calculations here, the modes are cut off 
abruptly, but an analysis of the shadow boundary transition would show a 
more gradual fading of each mode. For propagation from land to sea, the 


corresponding results hold. 


3) For propagation from sea to land, the land-tvpe groundwave modes excited 
at the shoreline are not detectable by an elevated observer who is too 
close to the shoreline. This shadow zone extends farther over land from 
the shoreline the higher the observer and the lower the order of the 
groundwave mode. For purposes of calculations here, the shadow zone is 
considered to have a sharp boundary for each groundwave mode, but an 
analysis of the shadow boundary transition would show a more gradual 


change. For propagation from land to sea, the corresponding results hold. 


4) Close to the shoreline, there is a direct edge-diffracted wave that con- 
tributes to the total field. For sea-to-land propagation, the edge- 
diffracted wave is an alternate representation of the field to the 
land-tvpe groundwave, and is more accurate near the shoreline. For 
land-to-sea propagation, the corresponding results hold. Which of these 
waves dominates the received field strength depends on the height of the 
observer and his distance from the shoreline, as shown on the accompanying 


plots. 


The results presented here apply not only to a land-sea boundary, but to 


any abrupt change in the surface impedance of the ground. 


3. THE CREEPING-RAY REPRESENTATION OF GROUNDWAVE MODES 


Keller (1962) showed how surface-diffracted waves could be represented by 
creeping rays (rays that follow geodesics on the surface). Levy and Keller 
(1959) derived the diffraction coefficients necessary to calculate field 
strengths for surface-diffracted waves. Jones (1968) recognized that groundwave 
modes on the Earth correspond to surface-diffracted waves and showed the equiva- 
lence of the surface-diffracted waves with the groundwave-mode sums of Watson 
(1918), Bremmer (1949), and Wait (1960). He also showed that representing the 
creeping rays as rays that have a particular complex height (that depends on 
the groundwave-mode propagation constant) gives a closer representation of 
the groundwave modes than having the creeping rays on the ground. Table 1 
summarizes the results of Jones (1968) and defines all of the terms it uses. 

His groundwave excitation coefficient and shedding coefficient together replace, 


but are equivalent to, the diffraction coefficient of Levy and Keller (1959). 


Table 1. 


Propagation along the ground 


Each groundwave mode has a cha- 
racteristic angular propagation 


a] a 
constant v. (e e 


gives the change 
in field strength of the ground- 
wave mode in a central earth an- 
gle a, and Vis complex to give 

both amplitude and phase.) It can 
be represented by a surface dif- 
fracted ray at a complex height 


v/k-a above the ground. 


Excitation of the groundwave 


The groundwave is excited by a 
ray tangent to the surface dif - 
fracted ray. The excitation coef- 
ficient equals the ratio of the 
transverse electric field of the 
groundwave at the point of ex- 
citation to that of the radiation 
field of the source at the same 
point. is the distance from the 
source to the excitation point. h 
is the height of the source above 


the ground. 


Radiation from the groundwave 


The groundwave sheds ravs (ra- 
diates) tangentiallv from the sur- 
face diffracted rav. The shedding 
coefficient equals the ratio of the 
transverse electric field a dis- 
tance £ from the point of shedding 
to that of the groundwave at the 
point of shedding exclusive of the 
phase integral and azimuthal 
focusing contributions (that is, it 
includes only elevation focusing), 
h is the height of the observation 


point above the ground, 


Propagation, excitation, and radiation properties of the groundwave 


Angular propagation constant 


VS ka «(E yi t 

where 

k is the propagation constant of free space. 
a is the radius of the earth. 

t is a root of Wi! (t) = q Wy (t) 

w(t) - Yr (Bi (t)- i Ai (t)) 


Ai and Bi are Airy functions 


TE YEN 


c is the surface impedance of the ground. 


is the impedance of free space. 


is the (complex) propagation constant of the ground. 


Excitation coefficient 


ij; 
E w,(0)/w,(1)? eit/4 (-y)!/4 12/3 (-y) / 


Dite- 2 w (y) 
t-q 
where ty i = kt TET rP, s! 
Special cases: 
Dy pr y(t) die w, (0)/w, (1) 
for small £, 
å : : : -lk 
Ditta Se w4(0)/w, (8) de ja EYE pikt,-ivtan vew (2) 1/34) 
t-q 
Z 2 
where h = (x) HA sas poat E 4? 
jo oye- 2 w, coy/w, y? (52)! /6(5)1/2. cin/4 ,, (o) 
-q 


t 


Shedding coefficient 
$ 3/2 
S(t) «(JE (ay l4 2/3 CY "y (yy fw, (0) 


where $65 = kt an t4 


special cases: 


s(£) ue us DRE s us jw, (0) 


lim 
t-00 


T S siek 
; e S)» elk f ivtan 77 tu t-(Gå) kn) jw, (0) 


j ES V 42 2 ey 
where h= (k) + £ od ot 


S(0) = 1 


After Jones, 1968; corrected 8 May 1972. 
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Physically, a groundwave mode is a ground-reflected wave for which the 
angle of incidence is such that the reflection coefficient is infinite. An in- 
finite reflection coefficient allows the boundary conditions to be satisfied by 
a reflected wave in the absence of an incident wave; therefore, a groundwave 
mode can propagate away from where it was produced. It is a reflected wave that 


satisfies the boundary conditions without an incident wave. 


Mathematically, groundwave modes arise from a residue expansion of a con- 
tour integral of a ground-reflected wave. The angle of incidence is related 
to the independent variable of integration in that contour integral; the poles 
of the integrand for the residue expansion are the poles of the ground 
reflection coefficient. The infinite ground reflection coefficient connects 


the mathematical with the physical interpretation of a groundwave mode. 


Figure 3 shows the geometry for calculating excitation of and radiation 
(shedding) from a groundwave mode. Jones (1968) shows that the vertical electric 


field at the receiver is 


-ik£ 


E ES T T Es cosB D (2 ) "DON ERI © 
4T an T vo 
XN sing 
—ik£ Ek T 
S(2) e v 55 osB, R (3.1) 


where the 7 factors are, respectively, (1) the incident field from the transmitter 
at the point of excitation of the groundwave, (2) the groundwave excitation 
coefficient D ps (3) the phase integral of the surface-diffracted ray from the 
point of excitation to the point of shedding, (4) the shedding coefficient S(2), 
(5) the phase integral from the point of shedding to the receiver, (6) a conver- 
gence factor due to azimuthal focusing, and (7) the pattern factor of the 
receiving antenna (a vertical dipole). The excitation and shedding coefficients 
are given in Table 1. Figure 3 can be used to evaluate all of the geometrical 
quantities in (3.1) in terms of known quantities. The quantity Ids is the current 


moment of the vertical dipole transmitting antenna. 


In the following sections I will make similar calculations to calculate mode 


coupling at the land-sea boundary. 


Receiver 


Center of Earth 


Surface of Earth 


Transmitter 


Figure 3. Representation of a groundwave mode by a surface-dif f'raeted ray at 


a complex hetght v/k-a. The sphere of complex radius v/k ts a caustic for 
the groundwave mode. From Jones (1968; Fide. 3). 


4. COUPLING OF GROUNDWAVE MODES BY EDGE DIFFRACTION 


As was mentioned in the introduction, over the sea we have one set of 
groundwave modes and over the land another set. A groundwave mode incident 
from the sea on the sea-land interface will excite groundwave modes appropriate 
to propagation over the land. One possible mechanism for conversion of a sea- 
type groundwave mode to a land-type groundwave mode is edge-diffraction. As we 
shall see, the hypothesis that edge-diffraction is the correct mechanism leads 


to agreement with independent calculations. 


Briefly, the edge-diffraction mechanism for coupling groundwave modes 
works as follows. The field of the sea-type groundwave mode incident on the 
land-sea interface is calculated in the usual way as though there were a re- 
ceiver there, but ignoring the pattern factor of the receiving antenna. The 
field strength is then assumed to be that of an incident wave on the land-sea 
interface at the appropriate angle of incidence. Using the appropriate 
scattering coefficient, we then calculate the scattered field from the land- 
sea interface. (This will be a cylindrical wave with a pattern factor 
determined bv the scattering coefficient.) We then treat that cvlindrical wave 


as a virtual source that excites a land-tvpe groundwave mode in the usual wav. 


Figure 4 shows a drawing of the ray structure for this process. However, 
Fig. 4 is inaccurate in that all four straight-line rav segments should be tan- 
gent to the circular surface-diffracted ravs. This is impossible to draw for 
the two rays that connect at the land-sea interface, even though the configura- 
tion is possible with complex rays. Figure 5 shows another representation of 
this process. Although it represents the surface-diffracted waves by creeping 
rays below the Earth's surface, and seems to require rays to penetrate the 
surface of the Earth, it correctly shows the four straight-line rays tangent 
to the creeping rays, and it can be used to calculate coupling of groundwave 


modes at the land-sea interface. 


We again consider a vertical electric dipole source. Then following the 
example in Section 3, the vertical electric field as seen by the observer in 


Ts. 15 


Free space 


Eo: Ho 


Observer 
Source 


Figure 4. Geometry to represent coupling of groundwave modes at a land-sea 
boundary by edge-diffraction. The impossibility of drawing a ray from 
the Land-sea interface that is tangent to the caustic for the groundwave 
mode disqualifies this figure as being useful for making calculations. 
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Free space 


Eo: Ho 


Z1 
£2 


Sea (^ Land 


E1 HOT 
TM 10:72 


Observer 


Source 


Figure 5. Geometry to represent coupling of groundwave modes at a land-sea 
boundary by edge-diffraetion. Although the rays are unrealistically 
represented as being below the Earth's surface, the correct geometrical 
relationships allow this figure to be useful for making calculations. 
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Idsuou. T -iv40, 
E = |i cosB D 6) e Soft) 
AT z T Vo T s, 2 
-ikk -ik£ —iv.a 
e * D4(-8, 785) peo he 1 D, Gj)| [e s 
vkk, 1 
-ikt, 
S m e (4.1) 
Vy R 


The 14 factors in (4.1) are respectively 


1) The incident field from the vertical electric dipole source at the point 


of excitation of the sea-type groundwave. 
2) The groundwave excitation coefficient. 


3) The phase integral of the surface-diffracted rav from the point of ex- 


citation to the point of shedding. 
4) The shedding coefficient. 


5) The phase integral from the point of shedding to the land-sea inter- 


face. 


6) The edge-diffraction coefficient for an incident grazing angle -B5 and 


a scattering angle -B> 


7) The geometrical spreading factor for the cylindrical wave scattering 


from the land-sea interface. 


8) The phase integral from the land-sea interface to the point of excita- 


tion of the land-type groundwave mode. 


9) The groundwave excitation coefficient. 
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10) The phase integral of the surface-diffracted ray from the point of ex- 


citation to the point of shedding. 

fi) The shedding coefficient. 

12) The phase integral from the point of shedding to the observer. 

13) A convergence factor to account for azimuthal focusing. 

14) The pattern factor of the receiving antenna (a vertical dipole). 

Table 1 gives the formulas for the excitation and shedding coefficients; 
I is the peak current in the source antenna; ds is the length of the source 
antenna; and the formula for the edge diffraction coefficient D. is given in 
Section 5. All of the other quantities in (4.1) can be found from the geometry 


in Fig. 5 in terms of known quantities. These are given in Appendix A. 


We can use Table 1 and (A.1) through (A.10) from Appendix A in (4.1). 


After some algebra, this gives 


fee s A MEN -t, \1/4 rales) y-t, 1/4 
pv, OLG 7 Ler 


V ? 172 Vo 3/2 Idswn 


Hr y ib 2i mad 


-iv,0 —-iv,0 
x ud " lu Wy (t5 -y4) w (t4-y) 


2 ka w, (t4) 


e 


(4.2) 


t t 2) 


— 2 — 
yas 1H 


where Appendix B gives the formulas for Yy» Y» V5» V» Io» d4» t5; and tj- 


We can make some approximations because the first term in (B.3) and in (B.4) 
is much larger than the second term. [At HF (10 to 30 MHz) ka is about 10°, and 


t, and t1 vary from order unity to about 50 for the first two hundred groundwave 
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modes. For high enough mode number, the approximation would eventually not be 
valid, but when that happens, (4.2) will no longer be valid because its 
validity depends on the first term in (B.3) and (B.4) being much larger than 
the second term.] With this approximation, the second factor in Square 


brackets in (4.2) is approximately unity, and we can use 
Vi - Va ~ ka (4.3) 


nearly everywhere but in exponential terms. Thus, (4.2) is approximately 


E ats) | Int, ju: "Y" y-ti |- Ge 
ži 2 o a 2 2m 2 r b 
k b “Vo kr “Vy 


2 1/3 1/2 
(2) (27) Dr 
—iv,9 -iv 0 
A 8.2 " ll w, (t-y,) Wi (t, y) 
2 2 OW EE.) w(t.) me 
t5-q, t4744 83 Er 


It is usual to normalize E to the free-space field above a flat, perfect con- 


ductor. We can take this to be 


Idswu —ika0 
o e 


MORI umet a Lade 4.5 
Ro $1 AT að ( ) 


Dividing (4.4) bv (4.5) and using (B.3), (B.4), and (A.11) gives 
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1/3/ y-t 1/4 1/3/ y-t 1/4 3/2, «3/2 
7 fk E 2 ka 1 a a 
DBT ES ka 2 z) | E J b. 2 ; | (2) (5) 
b ma, kr T 


siT/À 2 1/3 1/2 


-ix,t —ix.t 
B. tl w, (t,-y,) w, (t,-y) 


——————— ————— , (4.6) 
2 w (t5) w (t4) 


FI 


where X» and x, are defined in Appendix C. 
I will make more approximations after calculating the edge-diffraction 


coefficient in Section 5. 


5.  EDGE-DIFFRACTION COEFFICIENT FOR A LAND-SEA INTERFACE 


Figure 6 shows the geometry for edge diffraction by a plane interface in 
which the surface impedance changes discontinuously as it does at a land-sea 
interface. Let ur be the incident field at the land-sea interface that would 
exist if there were free space everywhere. Let u, be the diffracted field a 


d 
distance % from the land-sea interface. Let 


sinB + Å 
2 på 
R. ans (5.1) 
2 sinf, A, 
be the Fresnel reflection coefficient for reflection of the incident wave at an 


angle of incidence 1/2-B, from Medium 2. Let 


sinp- + A4 


R, = a m a (5.2) 
1 sing, A, 
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Free space no 


Uq 
4 
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NY 
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Medium 2 Medium 1 


Figure 6. Geometry for edge diffraction at a land-sea interfaee. Each of 
the two Lower media is characterized by a surface impedance 2. or z.. 


We consider a plane wave to be incident on the Land-sea interface. One 
ray of the cylindrical edge-diffracted wave ts shown 
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be the Fresnel reflection coefficient for reflection of a wave at the angle B1 
from Medium 1. A. Å,» B4» and B. are defined in Fig. 6. 


The incident wave u; gives rise to a diffracted wave 


e ikt 
u.D,(8,:8,) V j 


in the direction B. and a diffracted wave 


e ikt 
dg Ag) $ 


in the direction —BL- The latter gets reflected into the B1 direction to give 


a diffracted wave 


p ikt 


The incident wave also gets reflected to produce a wave 


u,R, « 


from the -B3 direction. This reflected wave in turn generates two diffracted 


waves, 
-ikk 
u. R Dp (B4 , 785) AC 
and 
-ik£ 
u. RD, (78, 3783) He > 


in the B1 direction. Thus, there are four contributions to the diffracted wave 


u Adding them together gives 


d' 


iT 


a ikt 
+ Zo. it Po) (5.3) 


Kaminetzky and Keller (1972) calculated the diffracted field for the geometry of 
Fig. 6. Their result [from their Equations (1.6) and (1.7) specialized to the 


case j=l] in the present notation is 


u. (2/0) 1/2 e 17/4 sinB, sinB, (AL-A.) -ik2 
a: a 1 ed NA a e (5.4) 
d (sinB,-A,) (sinB,-4,) (cosB,-cosB,) E. oi * i 
Using (5.1) and (5.2), we can write (5.4) as 
ġ 1/2 2345/4 ED CLERO Å-A Likt 
uq = u. (2/1) e — 7 post; ud. veg (5.5) 
2 1 vike 


Comparing (5.3) with (5.5), we see that they are the same if we take 


(5.6) 


We see that D, is even in both of its arguments. Thus, we can write (5.3) as 


E 
p ikt 
u, u, (l4-R)(l-R) D (8,,8,) —— . (552) 
Mo s" 3 X REIS x 


The above development requires a little comment. First, I altered (5.4) 
slightly from the original formula given by Kaminetzky and Keller. Their result 
was valid only for 
(5,8) 


La | «« |A | $ 
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and thev made no distinction between Al and A, in the denominator of (5.4). With- 
in the limitations of (5.8), (5.4) is equivalent to the expression given by 
Kaminetzky and Keller. However, the agreement of (5.4) with the formulas of Hill 
and Wait (1981), which we shall see later, indicates that (5.4) may be valid 


without the restriction (5.8). 


Second, I included the reflection coefficients R, and R, in (5.7) but did 
not include them in the 6th factor in (4.1), because in (4.1) the wave incident 
on the land-sea interface was a surface-diffracted wave (i.e., a creeping ray or 
groundwave). A groundwave already has a reflection coefficient included in it. 
In fact, the expansion in a groundwave mode series is a residue expansion at the 
poles of the ground reflection coefficient. A ground reflection coefficient is 
already included in both the groundwave incident on the land-sea interface and 


the groundwave excited at the interface. Because groundwaves occur where the 


reflection coefficient has a pole, we have 

[R | EE: (5.9) 
and 

IR, | 531 (5.10) 


for groundwaves, so that (5.7) becomes 


ej 
um u,R,RoD, (8, , 85) V 4 (5.11) 


where RR, is already included in the groundwaves. 

We can consider some other forms for (5.6). If the wave incident on the 
land-sea interface is a groundwave, and if the edge-diffracted wave is a ground- 
wave, then (A.6) and (A.7) give the incident and diffracted angles. Substituting 
(A.6) and (A.7) into (5.6) gives 
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ded (es) ar 11/4 A,-A, | ii 
E 2 (2m)1/2 9475 ; å 


Using (B.3) through (B.6) in (5.12) gives 


1/3 -it/4 q-q 
p. (us) Bollo e. (5.13) 


Another special case arises if the incident wave on the land-sea interface is a 
groundwave but the edge-diffracted ray is scattered at some arbitrary angle Bj: 


Then we Can use (A.6), (B.3), (B.5), and (8.6) In (5.6) to get 


1/3 e i7/4 4-4, 


D = i tka/2) WII 


: (5.14) 


(21) efr 9 


= 2 (ka/2)?/ 3 (cosg-1) " (5.15) 


rr 
H 


and I have written the scattering angle as B instead of SE 


6. COMPARISON WITH THE FORMULAS OF HILL AND WAIT 


Hill and Wait (1981) calculate groundwave-mode coupling at the land-sea 
interface. It is appropriate, therefore, to compare (4.6) with their formula. 


The appropriate edge-diffraction coefficient D_ is given by (5.13), because it 


E 
corresponds to both the incident and edge-diffracted waves being groundwave 


modes. Substituting (5.13) into (4.6) gives 
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1/3 .y,-t, \1/4 1/3/ y-t 1/4 
E/E = raft) PR IES P Bigs A | 
o | 2 FÆLE. 2 ONE 


1/3 


2 
i sinO t,-t, 
-ix,t -ix,t 
3 2 7 ji i11 w (t-y,) w (tj-y) 
2 2 w, (t4) v (E ) 
t5-d» t1744 TU 174 


Equation (6.1) is valid for arbitrary source and observer heights. 


(6.1) 


If the source 


and observer are close to the surface of the Earth, such that the angles zn and 


BR (Fig. 5) are small, then 


(B.1) is approximately 


a 


1/3 I/3 
s fa Y, K (2) (ka + (22) t,- kb) = t, -(2/ka) yn, , 


(B.2) is approximately 
1/3 1/3 
MA LTD da ka e d E 1/3 
ti $a (2) (ka + e t4 kr) = ti (2/ka) kh, 


and the two factors in square brackets in (6.1) are both unity. 


With these approximations, (6.1) becomes 


el 


(6.2) 


(6.3) 


(6.4) 


—ix,t —ix.t 
għ HWEREX 474» * 2-2 å pus 
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where, from (6.3) 
y, = (2/ka) kn, , (6.6) 
and from (6.4) 
y = (2/ka)! kn. (6.7) 


One further approximation can often be made if the horizontal distance be- 


tween the source and observer are small enough that 
seo. B. (6.8) 


Then (6.5) becomes 


1/2 (4474) exp (-ix,t,-ix,t,) Wi (t,-y,) w (t,-y) 


qw vua A (6.9) 
BA w, (t5) w (t,) 


E/E, = (mx/i) 5 
(t.-t,)(t,-4, ) (t174; 


where x is defined in (C.3). Of course, (6.9) gives the effect of only one ground- 
wave mode from the sea exciting onlv one land-tvpe groundwave mode at the land-sea 
interface. The total field is found by summing over all sea-type and land-tvpe 
groundwave modes. 

exp (-ix,t,-ix,t,) wy eo? w,(t,-y) 


3 XR E. EX , (6:10) 


necp 
E JE 9 (Tx/i) (q,-q.) ) ) 
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Equation (6.10) agrees with the corresponding formula of Hill and Wait (1981; 
bg. 9). 
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The agreement indicates that edge-diffraction is probablv the mechanism for 


coupling of groundwave modes at a land-sea interface. 


7. THE DIRECT EDGE-DIFFRACTED WAVE NEAR THE LAND-SEA INTERFACE 


Now that we know that edge-diffraction is the mechanism for coupling ground- 
waves at a land-sea interface, we see that groundwaves incident on such a land-sea 
interface will excite edge-diffracted waves in all directions, and that it might 
be possible to detect such waves near the shoreline. Figure 7 shows the geometrv 
that allows the calculation of these edge-diffracted waves. It is similar to 
Fig. 5 in that the creeping ray (the caustic for the surface-diffracted ray) is 
shown unrealistically below the Earth, so that the ray from the caustic to the 


land-sea interface can be shown tangent to the caustic. 


Following the example in Section 4, we can calculate the vertical electric 
field, seen by the observer in Fig. 7, due to a vertical electric dipole source. 


This is 


-k sinB,, 


get ren [sing]. (Jv) 


—-ik£ 
2 1 -ik£ 
eeo Ef] 


The 11 factors in (7.1) are respectively 


1) The incident field from the vertical electric dipole source at the point 


of excitation of the sea-type groundwave. 
2) The groundwave excitation coefficient. 


3) The phase integral of the surface-diffracted ray from the point of ex- 


citation to the point of shedding. 
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Figure 7. Geometry to represent excitation of an edge-diffracted wave by a 
groundwave mode tnetdent on a land-sea boundary. One ray of the 
cylindrical edge-dtffracted wave ts shown. 
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4) 


5) 


6) 


7 


8) 


9) 


10) 


LE) 


The shedding coefficient. 


The phase integral from the point of shedding to the land-sea inter- 


face. 


The edge-diffraction coefficient for an incident grazing angle -B, and 


a scattering angle B. 
One plus the ground reflection coefficient for the land. 


The geometrical spreading factor for the cylindrical wave scattering 


from the land-sea interface. 
The phase integral from the land-sea interface to the observer. 
A covergence factor to account for azimuthal focusing. 


The pattern factor of the receiving antenna (a vertical dipole). 


Equation (7.1) can be written in terms of known quantities by using the 


geometrical relationships in Fig. 7 (given in Appendix A), the groundwave excita- 


tion and shedding coefficients in Table 1 (with auxiliary formulas given 


explicitly for the geometry of Fig. 7 in Appendix B), and the appropriate form for 


the edge-diffraction coefficient (5.14). The result (after some algebra) is 


be EN Int; M/4 M stc AMA tai ^ 1/2 
2 Sie a 2 227 L kb r 
2 2 


1/2 -iv,60 
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where I have used 


Wo (t,)/w,(t,) - 44 


ER, «CO 1 ES (7.3) 
1 vi (t,)/w, (t4) 4 
for the ground reflection coefficient, and the Wronskian 
' m ! = =i 
w (t)w, (t) we (t)w, (t) T1. (7.4) 


As in Section 4, when (4.3) holds, the second factor in square brackets is 
unity. In addition, if we use (4.3) nearly everywhere but in exponential terms, 


and divide by (4.5), we get 


1/2 
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As in Section 6, if the source and observer are close enough to the ground (within 
about 100 km) then the factor in square brackets is approximately unity. In 


addition, we can use (6.2) in (7.5) to give 
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Further, if the source-observer distance is not too large, we can use (6.8) in 


(7.6) to give 


-ix,t -ikL+ikao 
2.2 w, (t,-y,) 1 


17:5 94 3957. e 
177 


teff 12 tod w (t3) (ki) 


E/E, = $(kad) sing 


ed Bd (7.7) 
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The formulas presented here are not valid near the horizon because of sur- 
face diffraction effects (excitation of groundwaves). I tried to extend the 
validity of the results by writing the field as a contour integral that I could 
evaluate by a saddlepoint approximation to get the formulas presented here. 

Near the horizon, however, the saddlepoint was near a pole and near an endpoint 
in the contour integration. Because of that and other difficulties, I was not 
able to extend the validity of these formulas closer to the horizon. Although a 
different contour integral representation might work, it is not of great 
importance to do this because the groundvave mode representation works well near 


the horizon. 


8. DIRECT RECEPTION OF A SEA-TYPE GROUNDWAVE OVER THE LAND 


Figure 8 shows how it is possible to receive directly a sea-type ground- 
wave over the land. If the Fender is high enough, then the possible pille 
effect from the land-sea interface will be negligible. The point is that the 
presence of the land in Fig. 8 should not be significant because the ray does 
not touch the land. The signal seen by the observer should be nearly the same 


as that for the case of a homogeneous sea. 


The condition that the observer is close enough to the land-sea interface — 


so that the situation in Fig. 8 applies is 


-1 2 
Br. cos S cos ai^ (8.1) 
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Geometry showing how a sea-type groundwave mode can be observed over 


Figure 8. 
land by an elevated observer. 


In Fig. 8, the radius V,/k of the caustic associated with the surface-diffracted 
wave is represented as being real. Normallv, however, Vo is complex, so that the 
radius Vo/k is also complex and the angles in (8.1) are also complex. Thus, in 


general, it may be appropriate to substitute 


y y 
wl 2 1 ^2 
01 « real (o 2 cos i) (8.2) 


for (8.1). I admit that this is an ad hoc procedure. However, it seems to have 
about the right properties and can probably be justified with more careful 
analysis. If the angles in (8.2) are small (as they usually are) then (8.2) may 


be approximated by 


X, < real (vy-t, - /-t,)s (8.3) 


where x is defined in (C.2), y is defined in (6.7), and t, is related to Vo by 


(B.3) and is defined as one of the roots oF 07). 


In practice, as 91 (or x4) increases, the higher order groundwave modes will 
cease to satisfy (8.3) first and be cut off. As 04 (or x4) continues to increase, 
more and more of the high order modes will be cut off until only the lowest order 
mode exists, and then finally it too will be cut off as z (or x4) continues to 


increase. 
The higher the observer (the larger y is) the larger 04 (or x) will be for 
cut off of a given groundwave mode. When the observer is on the ground, all 


groundwave modes are cut off beyond the land-sea interface. 


Figure 9 shows the cut-off situation for which (B.2) or (8.3) is mor 


satisfied. 
9. THE SHADOW ZONE FOR A GROUNDWAVE MODE EXCITED AT A LAND-SEA INTERFACE 


Figure 10 shows an observer over the land receiving a land-type groundwave 


mode. As pointed out in Section 3, a groundwave mode is a ground-reflected wave 
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Figure 9. Geometry showing cut off of a sea-type groundwave mode for an elevated 
observer who is too far from the shoreline. 
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Figure 10. Geometry showing how the ray that represents a land-type groundwave 
mode propagates in a stratght line from the Earth's surface to the observer. 
The ray must intersect the Earth's surface on the land side of the land-sea 
interface to satisfy the boundary condition (an infinite ground reflection 


coefficient) for a Land-type groundwave mode. 


al 


for which the reflection coefficient is infinite. (That allows a reflected wave 
to exist with no incident wave except at the point of excitation.) For the 
existence of the land-tvpe groundwave mode in Fig. 10, it is necessarv that the 


ray representing it come from the land and not from the sea; i.e., 


y y 
-1 -1 1 
94 ^ cos m COS va (9.1) 


Because the quantities in (9.1) are complex, we must substitute 


y V 
-1 1 -1 
0. > real [cos ES cos Ea (9.2) 


for (9.1). If the angles in (9.2) are small (as they usually are), then (9.2) is 


approximately 


X 


1 > real (vy-t, - Y-t} , (9.3) 
where X, is defined in (C.2), y is the normalized observer height defined in (6.7), 


and t1 is related to MI by (B.4) and is defined as one of the roots of (B.8). 


Figure 11 shows the opposite case, where (9.2) and (9.3) do not hold. The 
land-type groundwave mode seen by the observer in Fig. 11 appears to come from the 
sea. It is clear that the physical conditions necessary (an infinite reflection 
coefficient) for the existence of that groundwave mode do not hold there. It 
therefore seems likely that (9.3) is approximately a correct condition for the 
observation of the groundwave mode in question. There thus seems to be a shadow 
zone for each land-type groundwave mode. The observer in Fig. 11 is in the 
shadow zone for the groundwave mode represented, whereas the observer in Fig. 10 


is not in the shadow zone. 


In practice, an observer directly above the land-sea interface is in the 
shadow zone for all of the land-type groundwave modes that have been excited at 
the land-sea interface. As the observer moves over land, he moves out of the 
shadow zone of the higher order modes. If he moves far enough away from the sea, 


he will move out of the shadow zone of all of the land-type groundwave modes. 
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Figure 11. Geometry for an observer in the shadow zone of the land-type ground- 
wave mode that was exetted at the Land-sea interface. The ray that would 
connect the caustie of the Land-type groundwave mode with the observer 
tntersects the Earth's surface over the sea. At that point of interseotion, 


the raj fatls to satisfy the boundary condition (an infinite ground reflec- 
tion coefficient). 
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The closer he is to the ground, the sooner he will move out of the shadow zone of 
each mode. An observer on the Earth's surface is out of the shadow zone of all 


land-tvpe groundwave modes if he is over land. 


Although the analvsis presented here is somewhat ad hoc, it is probably a 


good approximation. 


10. THE FIELD NEAR A LAND-SEA INTERFACE 


There are three contributions to the total field observed near a land-sea 


boundarv. 


l3 Edge-diffracted rays from the land-sea interface (discussed in Section 


y P 
2) Sea-type groundwave modes (discussed in Section 8). 
3)  Land-type groundwave modes (discussed in Section 9). 


As pointed out in Section 8, the sea-type groundwave modes can still be 
observed for a distance over the land away from the land-sea interface. As 
pointed out in Section 9, the land-type groundwave modes excited at the land-sea 
interface cannot be observed by an elevated observer if he is too close to the 
land-sea interface. Thus, the significant variation of the groundwave mode 
contributions to the field is determined by the number of contributing modes 


for each of the two types. 


Figure 12 shows a plot of these three contributions for an observer 30 m 
above the ground [relative to the free-space field above a flat, perfect con- 
ductor (4.5)]. The edge-diffracted wave is very small near the land-sea 
interface, and increases as the observer moves away from the sea. The main 


cause of the rapid increase is the factor t in the denominator of (7.5). 


ag 5 
eff 2 
t, is small and complex with a postiive real part. As can be seen from (5.15), 


"abé 


is real and negative. The denominator taff t? will be smallest for small 
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Figure 12. The field (relative to the free-space field above a flat perfeet 
conductor) near the land-sea interface for an observer 30 m above sea level. 
The source ts a vertical electric dipole on the surface of the sea 20 km 
from the shoreline. The frequeney is 30 MHz. The numbers near the sea-type 
groundwave-mode contribution indicate the maximum mode number ineluded in 
the calculation. The numbers near the Land-type groundwave-mode contribu- 
tion indicate the maximum mode number not ineluded in the calculation. A 
maximum of 200 modes were used in these calculations. The dashed-line 
extension to the land-type groundwave-mode contribution indicates the sum 
of the first 200 modes. The dashed-line portion of the edge-diffracted 
wave indicates the region where the approximations used to caleulate the 
edge-diffracted wave are no longer valid because surface diffraetion by the 
Earth is significant. The values of to used tn (7.7) are given on the curve 


for the edge-dtffracted wave. Large negative values for t, (such as -50) 


indicate that the effect of surface diffraction by the Earth ts not 
significant. Because the edge-diffracted wave and the land-type groundwaves 
are alternate representations of the same field, the correct contribution 

to that part of the field will follow the edge-diffracted wave near the land- 
sea interface and the land-type groundwave far from the land-sea interface. 
The total field is found by summing (coherently) the sea-type groundwave 
with the edge-diffracted/land-type groundwave contribution. 
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B, that is, for edge-diffracted rays that are nearly horizontal. Thus, for a 
given height, that denominator will be smaller for larger distances from the 
interface, leading to the increase in signal strength shown in Fig. 12. The 
edge-diffracted rays are cut off at the horizon, but before that happens, the 
edge-diffracted ray solution is not valid because of surface diffraction (i.e., 
excitation of groundwaves becomes significant). The dashed part of the edge- 
diffracted ray contribution is probably not valid because of surface diffraction 


effects. 


The sea-type groundwave mode contribution in Fig. 12 continues essentially 
unchanged for about 500 m over the land. During that 500 m the higher order 
modes have been cut off because they no longer satisfied (8.3). At 500 m from 
the land-sea interface, only modes 1 through 15 have not been cut off. That 
the sea-type groundwave contribution has decreased very little during that 500 m 
indicates that only the first 15 modes were contributing significantly; the rest 
had been attenuated significantly during the propagation over the sea. During 
the next 500 m over land, the sea-type groundwave contribution falls rapidly as 
more and more modes are cut off. From 1000 m to 1200 m only one sea-tvpe ground- 
wave mode is not cut off. Beyond 1200 m all of the sea-type groundwave modes are 
cut off for an observer 30 m above the land. I suspect that a more thorough 
analysis would show the same qualitative behavior, but the distance from the 


shoreline where the rapid falloff occurs might be shifted a hundred meters or so. 


The land-type groundwave mode contribution in Fig. 12 is negligible at the 
shoreline and grows as the observer moves farther inland and out of the shadow 
zone of more and more modes. The observer moves out of the shadow zone of the 
higher order modes first. Finally, at 1200 m from the shoreline the observer 
moves out of the shadow zone of the lowest order mode, and beyond that distance 
all of the land-type groundwave modes contribute. The contribution of all of 
the modes calculated (200) is shown with a dashed line for reference. Again, I 


suspect that a more thorough analysis would show the same qualitative behavior. 


The agreement of the edge-diffracted ray contribution with the land-type 
groundwave mode contribution is significant. In fact, the two are not indepen- 
dent contributions, but rather alternate representations of the edge-diffracted 


contribution. The edge-diffracted ray is the more accurate representation 
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closer to the shoreline, while the groundwave mode representation is more 


accurate far from the shoreline. 


The practical importance of Fig. 12 is that an observer 30 m above sea 
level can be as much as 500 m back from the shoreline before his received 


Signal decreases significantly from what he would receive above the shoreline. 


Figures 13 through 19 show similar behavior for observer heights varying 
from 10 m through 300 m. Figure 20 shows similar behavior for a source 70 km 


from the shoreline. 


Figure 21 shows a contour plot in a vertical plane of the signal strength 
over the land for sea-to-land propagation. The three regions, the region 
where sea-type groundwaves are dominant, the transition region where some 
modes for both types of groundwaves are in the shadow, and the region where 
the land-type groundwaves are dominant are clearly shown. The exact shape of 
the contours in the transition region is uncertain because the exact transition 
of each groundwave mode across its shadow boundary was not calculated. The 
location of the transition region is reasonably accurate, however. Figures 


22 and 23 show more details of the contours closer to the shoreline. 
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Figure 13. The field near the land-sea tnterface for an observer 10 m above sea 
level. The other condtttons are the same as in Fig. 12. 
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Figure 14. The field near the land-sea interface for an observer 20 m above sea 
level. The other conditions are the same as in Fig. 12. 
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Level. 


The other conditions are the same as in Pigs 12, 
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Figure 16. The field near the Land-sea interface for an observer 50 m above sea 
Level. The other conditions are the same as in Fig. 12. 
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Figure 17. The field near the land-sea interface for an observer 100 m above sea 
Level. The other conditions are the same as in Fig. WIR, 
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Figure 18. The field near the Land-sea interface for an observer 200 m above sea 
Level. The other conditions are the same as in Fig. 12. 
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Figure 19. The field near the land-sea interface for an observer 300 m above sea 
level. The other conditions are the same as in Fig. 12. 
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Figure 20. The field near the Land-sea interface for an observer 30 m above sea 
Level. The source is a vertical electrie dipole on the surface of the sea 
70 km from the shoreline. The other conditions are the same as in Fig. 12. 
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Figure 21. Contours of constant field strength (relative to the free-space field 
above a flat, perfect conduetorjover Land near a land-sea interface. The 
source is a vertical electrie dipole on the surface of the sea 20 km from 
the shoreline. The frequency is 30 MHz. The conditions are the same as 
those for Figures 12 through 19. The dashed lines separate the three 
regions: the region where sea-type groundwaves are dominant, the transttton 
regton where some modes for both types of groundwaves are in the shadow, 
and the region where the land-type groundwaves are dominant. The exact 
shape of the contours in the transttton regton ts uncertain because the 
exaet transition of each groundwave mode across its shadow boundary was not 
calculated. 
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Figure 22. Contours of constant field strength over land near a land-sea inter- 
face, as in Fig. 21, but showing more detail near the shoreline. 
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Figure 23. Contours of constant field strength over land near a land-sea inter- 
face, as in Fig. 21, but showing more detail near the shoreline. 
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APPENDIX A. EQUATIONS FROM THE GEOMETRY OF FIGURES 5 AND 7 


2 gå 
hen b s (v,/k) ] (A.1) 
t = [a^ - ME 7 (A.2) 
t = få - EG (A.3) 
ka vm” (A.4) 
&, = cos (v,/kb) (A.5) 
8, = cos ^ (v,/ka) (A.6) 
B, = cog * (v4 /ka) (AZ) 
B, = cos (v,/kr) (A.8) 
04 = 04 — B1 - BR (A.9) 
a, = 0) - Ba - By (A.10) 
Ə = 0, +8, (A.11) 
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(APPENDIX A--Cont. 


) 


H^ +2 ar sin'9,/(1 É cos. )]^ 


P. m (a siné,/%) 


m/2 -o- 8 


1 


ejt 


/2 


(A.12) 


(A.13) 


(A.14) 


APPENDIX B. EQUATIONS FROM TABLE 1 APPLIED TO FIGURES 5 AND 7 


2 3/2 21.2 2 -1 

30% t,) = (kb - Vo ) - V, cos v,/kb, [B. 1) 

på ETEN bk ; Aa ee 

30 1 1 Vi cos 7 v,/kr, UB. 

/3 
å ka 

v F ka + 2 to» (BS) 
k 1/3 

Vi = ka + = zm (B.4) 


qom (B.6) 
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Z, and Zz, are the surface impedances of the sea and land, no is the impedance of 


free space, and t and t4 are determined by 


w,'(t,) = dow, (ty) Qu 
and 


w,'(t,) = quw, (ti) (B.8) 
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APPENDIX C. 
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DEFINITIONS 


(C.T) 


(C.2) 


(C.3) 


